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I. INTRODUCTION 



Muonium is one of the best studied purely electrodynamic bound states. The hyperfine 
splitting (HFS) in the ground state of muonium is measured with error bars in the 

ballpark of 16-51 Hz. A new higher accuracy measurement of muonium HFS is now planned 
at J- PARC, Japan [3]. The results of QED calculations of the HFS interval are usually 
organized in the form of a perturbation theory expansion in a, Za, m e /m M . Some of the 
terms in this expansion are enhanced by large logarithms of the fine structure constant 
and/or electron-muon mass ratio. The current theoretical uncertainty of the HFS interval 
is estimated to be about 70-100 Hz, respective relative error does not exceed 2.3 x 10~ 8 
(see discussions in (4 6 ] ) . Still unknown three- loop purely radiative corrections, three-loop 
radiative-recoil corrections, and nonlogarithmic recoil corrections (see detailed discussion 
in 5|, |6j) are the main sources of the theoretical uncertainty. Measurement of the HFS in 
muonium is currently the best way to determine the value of the electron-muon mass ratio. 
The value of a 2 (m^/m e ) is obtained from comparison of the HFS theory and experiment 
with the uncertainty that is dominated by the 2.3 x 10~ 8 relative uncertainty of the HFS 
theory jg]. Improvement of the HFS theory would allow further reduction of the uncertainty 
of the electron-muon mass ratio. A detailed analysis {4], 0] shows that reduction of the 
theoretical error of HFS theory in muonium to about 10 Hz is a realistic goal. As a step in 
this direction we consider below three-loop radiative-recoil contributions to HFS generated 
by the light-by-light (LBL) scattering diagrams in Fig. [T] (and by three more diagrams with 
the crossed photon lines). The radiative-recoil corrections due to the LBL diagrams in Fig. 
[1] are additionally enhanced by the large logarithm of the electron-muon mass ratio. The 
logarithm squared contribution was calculated long time ago 3]. Below we calculate the 
single- logarithmic radiative-recoil contribution. 
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FIG. 1. Diagrams with light-by-light scattering block 
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We will follow the general approach to calculation of three-loop radiative-recoil corrections 

on 

to HFS developed in |7H14I| and start with the general expression for the LBL scattering 
contribution in Fig. [TJ (see, e.g., j^, 5]) 



7T 3 F M\ 256 J J m 2 k A \k 2 + 2mk k 2 -2mk J W 11 



«7r 2 g 4 Vg 2 + 2Mq q 2 - 2Mg , 

where fc M is the four-momentum carried by the upper photon lines, g M is the four-momentum 
carried by the lower photon lines, m is the electron mass, M is the muon mass, Z — 1 is the 
muon charge in terms of the electron charged used for classification of different contributions, 
and S a pn U is the light-by-light scattering tensor. The Fermi energy is defined as 



= 3 M fe) '" C • < 2) 
where m r is the reduced mass. The angle brackets in Eq. ([I]) denote the projection of the 
7-matrix structures on the HFS interval (difference between the states with the total spin 
one and zero). 

The contributions to HFS of the first two diagrams coincide and with account for three 
more diagrams with crossed photon lines not shown explicitly in Fig. [1] we can represent 
the LBL block as a sum of two contributions, the first one corresponding to the first two 
(ladder) diagrams in Fig. [I] and the second one corresponding to the last (crossed) diagram 
in Fig. [1] 

f d 4 p 

S a /3au = / ~ n (2iwW + Cafjuu) , (3) 
J «7T 

where 



Lapp, = Tr ( 7 M - — lv — - — j p - — \ 7a — — j , (4) 

V ft ~ f ~ m v — m ft — f — m ft — m J 

\ ^ft—^—m ft — $ — ft — m u ft — ft — m a ft — m) 
The integral in Eq. ([1]) contains both nonrecoil and recoil corrections to HFS that are 
partially already calculated (see 0, [5] for a collection of these results) 



AE = ^£^[-0.472 514 (1)] + ™ ( £ In 2 — + C*! In — + C ) . (6) 
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The leading nonrecoil term in Eq. ([6]) is generated by the nonrelativistic pole in the muon 
propagator 

1 (7) 



and was calculated in 
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q 2 + 2Mq + iO M 



16| . This is a numerically dominant contribution and it should be 
extracted analytically from the expression in Eq. (pQ) before calculation of the radiative-recoil 
corrections. 

The diagrams in Fig. [T] contain three loop integrations and each of them could in principle 
generate a large logarithm of the electron-muon mass ratio. The strongly ordered region of 
integration momenta m<^ik<^ip<^iq<^iM would produce logarithm cubed contribution 
but it turns into zero due to the tensor structure of the LBL block and fermion factors in 
this region |s|, see below. The large logarithm squared, calculated in [?| arises from two 
integration regions, m<fc~p<Cg<M and m^k^p^q^ M. Our current goal is 
to calculate single- logarithmic contribution is Eq. (E]), and as a preliminary step we would 
like to separate the large logarithm squared contribution. First we recalculate the logarithm 
cubed and logarithm squared corrections in an intuitively transparent way. We will elucidate 
the origin of different contributions to the logarithm squared term what will help us to derive 
a convenient expression for calculation of the single-logarithmic contributions. 



II. LEADING LOGARITHMIC CONTRIBUTION 

We start calculation of the integral in Eq. ([!]) with the loop integration in the LBL 
scattering block. This block is a four-index tensor that is a function of four-momenta k 
and q. Tensor indices are contracted with the conserved antisymmetric electron and muon 
factors 



(8) 



As a result only the tensor structures that are odd in k and in q and are antisymmetric with 
respect to transposition of indices (a,f3) <H- (f3,a) and (//, v) -H- (z/, \i) give contributions to 
HFS. All tensor structures containing k a , hp, q^, and q u do not contribute to HFS. Then for 
our purposes the LBL scattering tensor depends only on two independent structures 

(9) 



•Map»v( k > <?) = \i.9^Kqp ~ 9np Kq a + g v pk^q a - g va k^qp). 



Due to symmetries of the electron and muon factors in Eq. (JB]) all terms in the tensor 
structures on the RHS in Eq. (Q give identical contributions to HFS, and in the calculations 
below we will use more compact expressions 

M^ik, q) = (k- q)g m g v $, M { ^(k, q) = g m k v q p . (10) 

The general expression for the LBL scattering tensor in Eq. fl3]) can be obtained by repre- 
senting the ladder and crossed traces in Eq. (J4j) and Eq. ([5]) in the form 

a ^ u D%D 2 D 3 ' a ^ u D^D^ 1 } 

where 



D 1 =p 2 -m 2 , D 2 = (p-q) 2 -m 2 , D 3 = (p - k) 2 - m 2 , D 4 = (p - q - k) 2 - m 2 . (12) 



Combining denominators in the expression for the crossed diagram with four different de- 
nominators in Eq. (1121) with the help of the Feynman parameters x, y, z, we obtain 



1 -x)D i + xl {l-y)D 3 + y 



zD 2 + (1 - z)D 1 



(p-K) 2 -n(x,y,z,0, (13) 



where 



K = k(l — xy) + q(l — x + xyz), 

Q(x, y, z,£) = m 2 — k 2 xy{\ — xy) — q 2 x{\ — yz)(l — x + xyz) (14) 
— 2(k ■ q)xy[l — x — z(l — xy)]£. 
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The additional parameter £ is equal one in Eq. ( TT3|) and is introduced here for further 
convenience, we will explain its role later. There is no denominator D± in the expression 
for the ladder diagram in Eq. (TIT]) but the expressions in Eq. (TT5T) and Eq. (TT4"]) at x = 1 
are still suitable for calculations with the ladder diagram. To obtain an explicit expression 
for the LBL scattering tensor in Eq. it remains to calculate the integral over the loop 
momentum and to write the result in terms of the independent tensor structures in Eq. (flOl) . 
The final expression for S a p^ v is rather cumbersome and will be presented below when we 
will use it for calculation of the single-logarithmic contribution in Eq. fl6]). We do not need 
that general expression for discussion and calculation of the logarithm squared terms below. 



A. The Region of Strongly Ordered Momenta 

Large logarithmic contributions to HFS arise from logarithmic integrations in the region 
of strongly ordered momenta m ^ k ^ q <^ M where the expression in Eq. ([1]) simplifies 

Arn a 2 (Z«) m n / 3\ f d 4 k h a h P ) f d 4 q M 2 q 2 . u „. n 

AE ~ \ T7 E F -77 / ■ 2 ,4 / ~^h~A 1777^ 7^7 )S a ^. 15 

7r 3 M V 64 / J nr 2 k 4 k 2 J %ix 2 q A q 4 - AM 2 q 2 ? ^ 

The factor with the large mass M is of order one when q <C M and integrations over q and 

k are logarithmic only if the LBL scattering tensor supplies a factor k/q. There are many 

ways how such factor arises in the expression for S a pn V . The leading contribution of this 

type could arise from the logarithmic integration over the loop momentum p in the region 

m<A;<p<g<M. (16) 

The integrand in the expression for the crossed diagram in Eq. (jSJ) contains large q 2 in 
the denominator in the region of strongly ordered loop momenta what makes integration 
over q nonlogarithmic. Only the ladder diagrams in Fig. [TJ could generate logarithm cubed 
contribution in this region. Expanding over q the ladder contribution to the LBL scattering 
tensor S a pp V in Eq. (j3]) we obtain 



= I — 2 L^ = [ d >Tr ( 7|t -— J lv-^l P —4 % 



■iTT 2 J in 2 \ ■ p-q-m f — m 'p — ^ — m p — m y 

d A p 1 ' ' 

-7^8 Tr (MlvtlPfttlo 



Projection of the last trace on the structures in Eq. (Q is zero after averaging over di- 
rections of vector p^. Hence, leading term in the expansion of the LBL tensor does not 
contain (k/q) ln(q/k), and the logarithm cubed contribution to HFS does not arise despite 
the logarithmic nature of all integrations in the region of strongly ordered loop momenta in 
Eq. (US]). 



B. Logarithm Squared Contribution 

The logarithm squared contributions arise in two integration regions 



m<A;~p<o-<g<M, m</c<cr<p~g<M, (18) 

that are obtained from the region of strongly ordered loop momenta in Eq. (flu]) , when we 
lift the strong ordering requirement and allow two of the three loop momenta to be of the 
same order. For calculational purposes we introduce an auxiliary parameter m <C a <C M 
that separates the regions of large and small momenta and should cancel in the final results. 
As we have seen the LBL scattering tensor does not generate leading logarithmic terms of 
the type (k/q) ln(q/k), and hence logarithms (k/q) ln(a/k) and (k/q) hx(q/a) do not arise in 
the large and small integration momenta regions in Eq. (fl8|) . 

To isolate the logarithm squared contributions we expand ladder and crossed contribu- 
tions to the LBL scattering tensor in Eq. fl3]) in the regions of small and large momenta in 
Eq. ( fl8l) and look for contributions of order k/q. The ladder contributions of this type arise 
from expansions both in the small and large momentum regions 

[ d A p ( 1 1 1 1\ 

C a ^\p <a]^J —Tr [1,-1,-1,— Jla- ) , (19) 



J m \ p- 4 p p p pj 



(20) 



The low momenta integral in Eq. f fl9|) is superficially linearly divergent at the upper limit. 
This is an artificial divergence introduced in the integral when have thrown away momentum 
p in comparison with q in one of the denominators. We will deal with this divergence below 
extracting contributions of order k/q from the integral in Eq. ( Tl9l) . Notice also that in the 
large momenta region the integral in Eq. (|20|) generates a contribution of the form k/q after 
expansion of the integrand up to the second order in k/p. 



The contribution of the small integration momenta region for the crossed diagram in 
Eq. (JSJ) is suppressed as a 2 /q 2 due to large q 2 in the denominator of the crossed diagram 
in Eq. §5§. The total leading logarithmic contribution to HFS connected with this diagram 
arises only from the large momenta region, where the terms of order k/q arise, like in 
Eq. fl20p . after expansion of the integrand up to the second order in k/p 

fd^p^ ( 1 1 1,1 1 1 1 1 1 1\ , , 

c '"- * J l^ Tr V^f^y^y + ^jzrf^jzrfy^) • < 21 > 

Now we are ready to calculate the logarithm squared contributions to HFS. We rational- 
ize and combine denominators in Eq. f fl9|) -Eq. (j2ip like in Eq. (TT2|) -Eq. f JT4|) . and extract the 
contributions proportional to k/q from the respective integrals. Let us illustrate necessary 
transformations considering as an example calculation of the small momenta ladder contri- 
bution in Eq. ( fl9l) . In this case we need only one Feynman variable y in Eq. ({TBI , other 
variables are fixed, x = 1, z = 0. We obtain 

c a ^\p<a}~^ I I dy — — 2 y Tr [ 7 ^ 7 ^^_^) 7Q|) ]. 

q J ™ Jo [( p -k(l-y)) +k 2 y(l-y)\ 

\p\<a L J 

(22) 

Next follows shift of the integration momentum p — > p + k(l — y) in this formally linearly 
ultraviolet divergent integral. As usual a finite surface term arises after such shift (see, e.g., 



17(|). Further transformations are pretty standard, we calculate the trace, make the Wick 
rotation, preserve only the contributions proportional to k/q, calculate the integrals, extract 
the terms proportional to the tensor structures in Eq. (1101) . and obtain a finite result 

4 

£ a piu,\p <cr] = -^[(k- q)g lia gp v - 1g m qpK ] . (23) 

Notice that the coefficient before the formally linearly divergent contribution turned into 
zero automatically. Calculation of the contributions in Eq. ( 1201) and Eq. (|2~TT) is no more 
difficult, and we obtain 

4 

ZaPnAP > a } = -^g^akuqp, (24) 
g 

C a ^u = -^Qm [(k ■ q) g u f3 - %Kqp\ . (25) 
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Finally, the total contribution to the LBL scattering tensor Sap^ of the type k/q arising 
in the region in Eq. ffl8|) is the sum of the contributions in Eq. fl23l) -Eq. fl25|) 



16 

[p < cr] + 2C a p^\p > a] + Ca^v ~ — [(k • q)g m gpu - Zg^aqph] . (26) 

Comparing the small momenta contribution to the LBL tensor 2L a ^^ v \p < a] (see Eq. ( 123]) ) 
with the final expression above we see that the small momenta contribution to S a p^ v is 
equal the large momenta contribution. Next we substitute the LBL tensor in Eq. ( 126|) in 
the expression for HFS in Eq. ( TT5l) and contract Lorentz indices 

(l a h P )(l"il U )g, a gup = -~ [ 2(k ■ q) + Mo] , (27) 

(rh^il^g^Kqp = \ [k 2 q 2 - (k ■ qf + k 2 q 2 + fc 2 g 2 - 2(k ■ q)k q ] . (28) 

Averaging over directions of k we obtain the ladder and crossed contributions to the LBL 
tensor 



q o q 



(29) 



(i a h")d^Y)c a ^ = (i a h p )(r<Hi v )j 2 g» a [(* • g)fc* - 3^] = - A 2 ^±M. (30) 

Then the total LBL scattering tensor in this regime has the form 

( 7 ^){j^Y)(2C a ^ + C a ^) = -32k 2 2q2 + q \ (31) 
and the logarithm squared contribution to HFS can be written as 

a\Za) m 3 f d*q M 2 (2q 2 + q 2 ) f d*k 

7T 3 M F 2J ivrV g 4 -4M 2 g 2 J 1 ' 

The logarithm squared contribution arises in the integration region m 2 <C k 2 <C q 2 <C M 2 , 

and we calculate it using, after the Wick rotation, four-dimensional spherical coordinates 

and preserving only the logarithmic contribution 
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AE 



a 2 (Za) m 



7T 1 



M' 



E h 



3\ 2 

"2 / 7T 



d9 sin 2 6 

J ni- 



Al 2 



dq 2 M 2 (2 + cos 2 



g 2 q 2 + AM 2 cos 2 6 ] 2 ~P~ 



9a 2 {Za) m _ . 2 M 
— -£/pln — . 

M m 



(33) 



7T" 



This logarithm squared contribution was obtained in Q]. 

Another way [18| to calculate this contribution is to notice that it is intimately connected 
with the two-loop renormalization of the axial vector current first calculated by Adler long 



time ago 



191 ] . Consider the leading recoil correction to HFS generated by the graphs with 



two-photon exchanges in Fig. |2j Respective contribution to HFS is given by the expression 
in Eq. ([I]) without the LBL scattering block. The antisymmetric electron-line and and 
muon-line spin factors in Eq. (jHJ) can be written in the form 



(TVA.) — > ^( 7 p75>(e), (7*7*7%) — > «"^<7A7te> 



(/*)■ 



(34) 




FIG. 2. Diagrams with two-photon exchanges 



Then the leading recoil contribution to HFS has the form 



AE 



■(7a7s) 



(c)- 



(35) 



3 Za in . \ f M 
4TM £F(7 75) W 1 q 
We see that the leading recoil contribution to HFS is determined by the matrix element of 

the electron axial current calculated at the characteristic virtuality q. The first radiative 

correction to this matrix element is of order a 2 ln(g 2 /m 2 ) and arises at two loops 19| (see 

Fig. E} 



fx -> fx 



4 \7T/ 



(36) 



Substituting this renormalization factor in the expression for the leading recoil correction 
in Eq. (|35|) we obtain the leading recoil correction accompanied by the logarithm squared 
contribution in Eq. fl33|) 



3Za m M 9a 2 (Za) m 



AE = —— — ——-Epln h , 

4 ii M m 4 



7T" 



^Epln 2 —. 

M m 



(37) 
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7a Ys TATs 



FIG. 3. Two-loop axial current renormalization 

We see that the logarithm squared contribution is an observable effect of the axial current 
renormalization. 

Let us clarify the connection between the calculation of this contribution based on the 
consideration of the LBL scattering tensor and the approach with the axial current renor- 
malization. The diagrams in Fig. [3] are naively linearly divergent and the result in Eq. ( 13T)j) 
implies a gauge invariant regularization. On the other hand contribution of the diagrams in 
Fig. [I] to HFS is ultraviolet finite. Hence, for the purpose of calculation of the logarithm 
squared contribution replacement of the two triangle diagrams in Fig. [3] by the three box 
diagrams in Fig. [I] can be considered as a gauge invariant regularization. Let us see how 
this regularization works. Consider first the integration region where k ~ p < o < q. In 
this region the lower fermion line in the LBL scattering box in the first two diagrams in Fig. 
Q] effectively shrinks to a point, and these diagrams turn into the axial current diagrams in 
Fig. [3j We have already calculated contribution to the LBL scattering tensor 2C a p^ u [p < a] 
generated in this region, see Eq. ([23]) . This contribution is exactly one half of the total 
contribution to the LBL scattering tensor in Eq. ( 1261) . The other half (see Eq. (1241) and 
Eq. ( 1251) ) arises in the integration region k < a < p ~ q, where both lower loops in all 
diagrams in Fig. [1] effectively shrink to a point. Therefore, contribution to the axial current 
renormalization generated in this region can be considered as a pure regularization effect, 
for example contribution of a heavy regularizing fermion that survives when the fermion 
mass goes to infinity. 



III. SINGLE-LOGARITHMIC CONTRIBUTION 



Calculation of the single-logarithmic contribution to HFS requires more accurate treat- 
ment of the LBL scattering tensor in Eq. fl3]). Calculating traces in Eq. fTTTj) we obtain 
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x [(p • g) + (p • Jfe) + (A: ■ q)} + 32^ Q [(Jfe ■ q)p v pP ~ (p • g)fcV ~ (p • *OpV] ( 38 ) 
+ 16(p • fc)(p ■ q)g m g vfi - 32jf P a k v q^, 

and 



= 8£>i^ a [-3&V + fcV + pV] - 8D 2 g fia k u p^ - 8D 3 g^ a p u q^ 

+ 16(Jfe • q)g tia p v pP - 160 • q)g m k v pP - 16(ife ■ q)g'" x p v q p + 16p ■ (A; + q)g m k v q$ . 

(39) 

Then after the shift in Eq. ( Fl3|) and calculation of the loop integral the LBL scattering tensor 
can be written in the form 



(40) 



where (see Eq. (Til 



and 



7T 2 2 



1 r ± 



10 / / dz^ 7 ^——^^[{k-q)g m g v p-g m k v qp] 



+ 



7o ^C 1 ,?/,*, 1 ) 
Jo Jo I ^l 1 ^ 



M) 



+ y 2 (l-^ 



fc 2 (l - y) 2 + q 2 y 2 z 2 



W(l,y,z,l) 



(k ■ q)g m g v p - 2g lia k 

uqp 



+ 16 J dy [ dz 



m 2 y 



+ y i (l-y)z(l-z) 



_ f 1 d ^k-q) 2 y\l-y) 2 z 2 i 
o ~~ | [Q(l,y,z, 1) J Q 2 (l,y,z,^) 

k 2 q 2 + (k ■ q) 2 ' 



n 2 (l,y,z,l) J 



g^agu/s 



(41) 



QCt/3/lV 



d A P rial 



7T 2 i 



1 r l 



+ 



= 8 / dx dy dzx 2 y 
Jo Jo Jo 

Am 2 — 2k ■ qxy[l — x — z(l — xy)\ 



(k ■ q)g^ a g v p - 3g m 



tt(x,y,z, 1) 



Q 2 (x,y,z, 1) 



g m k v qa }■ 



(42) 
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The additional parameter £ in Eq. (HI]) arises when we separate the ultraviolet divergence 
in the logarithmically divergent integral 



d A p 



P 



+ 



m 2 [p 2 - Q(l,y,z, l)f 
d 4 p„ 2 



d 4 p 



TP 



m 2 [p 2 -tt(l } y } z,0)} 3 
1 



[\p 2 -n(l,y,z,l)\ 3 (p 2 -n(i, y ,z,o)r 



(43) 



In 



A 2 



y,z,0) 2 



1 2(k-q)y{l-y)z 



As was explained above the LBL scattering tensor is contacted with the odd in k and in 
g tensor structures and the even ultraviolet divergent term in the square brackets does not 
contribute to HFS and can be thrown away. 

It is easy to check that the terms on the RHS in the first lines in Eq. (|4"T|) and Eq. ( 142|) 
reproduce the leading ladder and crossed contributions to the LBL scattering tensor in 
Eq. (1261) and generate the logarithm squared contribution to HFS in Eq. ( 1331) . 

At the next step we calculate the integral over momentum k of the upper photons in 
Eq. © 



T(q 2 ,q ) 



1 



d 4 k 



1 



1 



2 ./ m 2 k 4 \k 2 + 2mk k 2 - 2mk 
d 4 k (i a h P ) 



(44) 



-s, 



in 2 k A k 2 — 2mko 

where we used the symmetry of the integrand under simultaneous substitution k — >■ —k and 
q — > —q to get rid of the second term in the first brackets on the RHS in the first line. In 
terms of T(q 2 , q ) the contribution to HFS in Eq. ([T]) has the form 



AE 



a 2 (Za) m 



7T" 



M 



3M 2 
128 



d\ 

ivrV \q 2 + 2Mq ' q 2 - 2Mq 



+ 



T(q 2 ,q ). (45) 



We calculate T(q 2 , q ) combining denominators with the help of additional Feynman 
parameters u and t (the four- vector Q and the scalar A on the RHS depend on the parameters 
x,y,z,£,u, and t) 



u 



4 - t)k 2 + t{k 2 - 2mk ) 



u 



-xy(l - xy) 



(k-Q) 2 -A, 



(46) 
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where 



A = g -g 2 + 2bq + a 2 
1 — x 



r 2 + 



m 2 u 



xy(l - xy) 



b = — , 
9 



d = 

go 



z — 



1 — xy 
u(l — yz)(l — x + xyz) 



T = m(l - u)t, g = g -d 2 



y(l - xy) 
The four- vector Q has the form 



(47) 



= dq^ + Tp, (48) 
where = (r, 0) and to calculate the integral over k we shift the integration variable 



kf, — >kft + Qp = k,j, + dq^ + Tp. (49) 

Contractions of matrix structures in Eq. (j44p are calculated using relationships in Eq. ( 127|) 
and Eq. ( 128]) . Five different tensor structures arise in calculations, and after the shift of 
integration variable in Eq. (1491) they reduce to 



8 



(i a h P )( f fff)(k'q)g m g vP 



(2q 2 + q 2 ) - 8q r, 

| (\k 2 + q 2 dA (2q 2 + q 2 ) - 8q 2 Q r 2 - jj(5g 2 + q 2 )dq r, 



kfh P )kfff)9^Xqp — ► ,k 2 (2q 2 + q 2 ) + *(q 2 - q 2 )r 2 
(l a h P )(l»il u ){k.q) 2 g m g uP 



3 - v i ■ iu/ ■ 3 



o ^ A k 2 + q 2 d 2 )(2q 2 + ql)q 2 d - "-(8q 2 + q 2 )r 2 q 2 d 



- 8 -( 1 -k 2 + q 2 d 2 )(7q 2 + 2q 2 )q r-8qy 



3 V4 

(l a h $ )(l^YKk-q)g m k v q p 



^k 2 (2q 2 + q 2 ) + °-r 2 (q 2 -q 2 ) 



q 2 d + q T^j . 



(50) 



Now we are ready to obtain explicit integral representations for the ladder Ti(q 2 , q Q ) and 
crossed Tc(q 2 ,qo) diagram contributions to the function T(q 2 ,q ) 



T(q 2 ,q ) = 2T L (q 2 ,q )+T c (q 2 ,q ). 
14 



(51) 



The ladder contribution can be written as a sum of nine integrals 



Th{i \ go) = \r- [ dy [ dz f du f dt ^ 7i,i(y, z, u, t, q 2 , q ), 
<J Jo Jo Jo Jo i 

where (x — 1 in all formulae below and £ = 1 in all expressions except 71,6) 



T L ,i = yz{l-t){l-uf 



1 q 2 d 2 

A A 2 



( 2g 2 + q l) - i q2 + 2%2)r2 q °( 5q2 + q ^ Td 



A 2 



A 2 



r 3 (o n 2 Mn 2,} (l-2y) + 2yz(l-t)(l-u) 2 
Il,2 = ~{2q + q ) 



+ (!-*) 



l-y A 

u(l-u) y 2 z 2 (l-z)q 2 (l-t)u(l-u) 2 



A 2 



Tj ^ (l-2y) + 2^(l-t)(l- M ) 2 _ (1 _ - u) 



+ 2 



l-y A 2 
yV(l-^)g 2 (l-t) M (l- M ) 2 



(i-y) 5 



A 3 



A 2 



■(2g 2 + %W, 



77 



LA 



(l-2y) + 2yz(l-t){l-u? 



A 2 



_ + 2 .^ 2 (! - ^ 2 (! - *)«(! - ^) 2 



x 



71 



A 2 (l-y) 2 
(2g 2 + g 2 )r 2 + g (5g 2 + g 2 )rrf 



A 3 



m 2 (1 - i)(l — u) 1 



71 



L,6 



l-y A 2 

= 4 f dttyz 2 {l-t)u{l-u) 2 
Jo 

1 ^2^2 



(2g 2 + g 2 )rf + 3g r], 



x 



+ 



3 1 g 2 d 2n 



4 A 2 A 3 



1 1 g 2 4 
4 A| ~ "Af 



(2q 2 + q 2 )q%- T -j^(8q 2 + q 2 ) 



(7g 2 + 2g 2 )g r-M^ 



7i, 7 = 



yz(l — z) q 2 u{l — u) 



i-y 



A 2 
15 



(2g 2 + g 2 )rf + 3g r], 



77 



L,8 



2 yz{1 - Z \ l-t)u(l-uf 



i-v 



3 1 q 2 d 2 



1 1 

~4A 2 ~ + A 3 



4 A 2 A 3 
q 2 d 2 



(2q 2 + ql)q 2 d 
3g 3 r 3 



A 3 



T f , 9 = 4 ^ (1 Z \ l-t)u{l-uf 



i-y 



1 1 
4A~ 2 ' 



-l^-(2q 2 + ql)q 2 d 



r 2 d 1 1 



+ q 2 (q 2 - <? 2 )^- - 1 a 2 ^ 2 ^ + q ^ q ° T + %< " q2 ~ q ^ A 3 

The crossed diagram contribution can be written as a sum of three integrals 



(60) 



(61) 



Tc(q 2 ,qo) = ^r dx dy dz du f dt^Tc^x^, z,u,t,q 2 ,q ), (62) 
<J Jo Jo Jo Jo Jo j 

where (£ = 1 in all formulae below) 



l x(l-t)(l-^ 2 

'C,l — o 1 

2 I — xy 



(V + % 2 ) 



2 g 2 c/ 2 ]_ 3 g 2 r 2 g (5<? 2 + g^)rci 



A A 2 



A 2 



A 2 



Tc,2 - 
Tc'3 



x(l -t)(l -u) 



um 



1 — xy xy(l — xy) 
I x{l - t){\ - u) 2 



2q 2 + ql _ 4 (q 2 - q 2 )r 2 



A s 



A 3 



1 — xy 



o 2 d 2 o 2 r 2 d 2 

(V + g 2 )^r-% 2 -% 2 ) g 



A 3 



+ (2q 2 + q 2 ) q ^-<q 2 -q 2 f y " 



(63) 
(64) 

(65) 



A 2 ^ A 3 

All logarithmic contributions to HFS can be obtained from the large momentum expan- 
sion of T(q 2 , q ). The leading term in this expansion we already obtained in Eq. ( I3T1) 



T 



-32 



d 4 k q 2 + 2q 2 
Z7T 2 /c 4 q 2 



-32 



- q dJf_ q 2 + 2q 2 
2 k 2 a 2 



-32 



q 2 + 2q 2 -q 2 



In 



We are going to calculate the next terms in the large q expansion of T(q 2 , q ) 



T = -32 



q 2 + 2 ql ,_ -t 



ln 



2g 2 + q 2 q 2 + 2q 2 

+ «2" 



(66) 



(67) 

q* m z. q« ^ 

where K\ and k<i are the numerical coefficients to be calculated. The numerators of the last 
two terms should contain any two linear independent combinations of q 2 and q\. We have 
chosen 2g 2 + q^ and q 2 + 2q 2 . because the first one contains the same combination of momenta 
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TABLE I. Leading terms in the ladder diagram expansion 



7l,i 


16 Kl « 


3 „ 
T6 K 2i 


Tl,i 


ln^ + 8C(3)-8 


1 


71,2a 


-61n^- -9 





71,26 


31n^-24C(3)-^ + 35 





71,2c 


2 

31n^-24C(3)+33 





71,3 


96C(3) - 116 





71,4 


4-7T 2 8 

9 3 





71,5 








71,6 


-16C(3) + 19 





71,7 


-16C(3) + 20 





71,8 


-40C(3) + f 





71,9 


16C(3) - 19 





71 


m z 9 6 


1 



that arises in the leading term of the expansion in Eq. f l66|) . and the second structure does 
not generate logarithm of the mass ratio after integration over q. 

The integral representations in Eq. ( )53|) -Eq. (!6Tl . and in Eq. ( )63|) -Eq. (|65|) are rather 
cumbersome and we simplify them before integration over the Feynman parameters. It 
turns out that in calculations of the asymptotic expansion in Eq. ( IBTl) we can omit the term 
2bqo in the denominator A 

A = g(q 2 + 2bq + a 2 ) g(q 2 + a 2 ) = A. (68) 

To justify this simplification it is sufficient to notice that the term 2bqo in the denominator 
A arises when we combine the subleading term 2mk Q from the electron propagator with the 
subleading term 2k ■ q from the denominator Q in Eq. fj46|) . The denominator A depends 
only on the Lorentz invariant momentum squared q 2 what makes the calculations easier. 
Even after this simplification the integrals over the Feynman remain unwieldy, especially in 
the case of the crossed diagrams when they contain an extra Feynman parameter x. They 
become more manageable if we notice that the remnant of the scalar product k ■ q in the 
denominator Q survives not only as the term 26g m Eq. (168]) but also as the second term 
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in the factor g (see definitions in Eq. (14 7p . below we show explicitly only dependence of g, 
d, and A on one parameter £) 



9(0 = 9(0) ~ d 2 



u(l — yz)(l — x + xyz) 



1 — X 
1 — xy 



(69) 



- xy) 

The leading logarithmic term in Eq. (I6"7j) arises when we simply omit this second term, 
g(£) — y g(0)- To calculate the sub leading terms in the asymptotic expansion in Eq. (1671) we 
represent all terms with the denominator A(£ = 1) in Eq. (168]) in the form 



1 



1 



1 , ^q 2 d 2 (0 , x 
d£ ' ■ (70) 



o 



A(0 



A(£ = l) A(£ = 0) ' [A(£ = l) A(£ = 0)J A(£ = 0) 

It is easier to calculate separately the integrals with the two terms on the RHS than the 
integral with the denominator A on the LHS side. Unlike the expansions in Eq. (I24]) -Eq. ( 126]) 
the integrals with the denominator A substituted by the terms on the RHS are well suited for 
calculation of the subleading terms in the asymptotic expansions in the kinematics described 
in Eq. ( 1T8]) . After tedious calculations we obtained analytic expressions for the subleading 
terms in the asymptotic expansion of the integrals in Eq. (I53p -Eq. ( 1611) and in Eq. (I6"3"]) - 
Eq. (163]) . These terms are collected in Tables UJ and [Til an d the asymptotic expansions for 
the ladder and crossed diagrams have the form 



16 2q 2 + q 2 Q 



In 



57T 



77T 



16 q 2 + 2q 2 



(71) 



T, 



c 



64 2q 2 + q 2 
' 3 q 2 



In 



2C(3) 



32 q 2 + 2q 2 



Then the total ultraviolet asymptotic expansion in Eq. (J67J) acquires the form 



(72) 



T 



-32 



2? 2 + ql 
q 2 



-q 2 4 
ln^--C3 
m z 3 



8^^ 

"27 



37 
18 



64 q 2 + 2q 2 
3 q 2 ' 



(73) 



Our next task is to calculate the single-logarithmic contributions generated by the first 
two terms in the expansion in Eq. ( 173]) . We substitute the asymptotic expansions in Eq. ( 1711 
and Eq. (172]) in the expression for HFS in Eq. ( 145]) . One can prove that the last terms 
with the numerator q 2 + 2q 2 do not generate logarithmic contributions and all double- and 
single- logarithmic contributions are generated by the terms proportional 2q 2 + q 2 . The 
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leading logarithmic term in Eq. ( I73|) generates not only the logarithm squared contribution 
in Eq. ( 133]) but also additional single-logarithmic terms. Analytic calculation of single- 
ogarithmic terms is performed with the help of auxiliary integration formulae collected in 



201 ] . We obtain ladder and crossed diagram logarithmic contributions to HFS in the form 



2AE T 



a 2 (Za) m 



7T" 



M 



E f 



3, 2 M 
-In 2 — 

8 m 



7T 2 23\ M 

1 In — 

3 16 / m 



AE, 



c 



a 2 (Za) m 



7T" 



M' 



3, 2 M 
-In 2 h 

8 m 



-3C(3) 



17\ , M 
— In — 

2 m 



(74) 



(75) 



TABLE II. Leading terms in the crossed diagram expansion 



"64 Kli 



3 



Tc,i 

Tc,2 
Tc,3 



ln^-^ + l 
-| C (3) + ^ + l 



|C(3) 



ii 



l 

2 






Tc 



ln=%- 



2C(3) + 



And finally the total logarithmic radiative-recoil contribution to HFS generated by the 
gauge invariant set of of three-loop diagrams with the LBL insertions has the form 



AE 



a 2 (Za) m 



TV" 



M' 



Ei 



4 m 



2tt 2 91 \ , M 
1 In — 

3 8 m 



(76) 



IV. CONCLUSIONS 

Other single- logarithmic radiative-recoil contributions to HFS were calculated earlier 



a 3 m 



M 



AE = ^r^-E F [3C(3) - 6vr 2 In 2 + tt 2 - 8l In — . 



(77) 



7T 3 M m 
Combining this contribution with the result obtained above in Eq. f J76|) we obtain the total 

result for all known three-loop radiative-recoil single-logarithmic corrections to HFS 



a 3 m 



7T 3 M 



7T- 



27 



AE tot = -,—E F -6n z In 2 + — + — In — 



M 



(78) 
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As was explained in the Introduction, the current goal of the HFS theory in muonium is 
to reduce the theoretical uncertainty below 10 Hz. The result above is a step in this direc- 
tion. Work on calculation of the remaining three-loop single-logarithmic and nonlogarithmic 
contributions to HFS is now is progress. 
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